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Abstract 

We show that the field equations for cosmological perturbations in Newtonian 
gauge always have an adiabatic solution, for which a quantity TZ is non-zero 
and constant in all eras in the limit of large wavelength, so that it can be 
used to connect observed cosmological fiuctuations in this mode with those at 
very early times. There is also a second adiabatic mode, for which TZ vanishes 
for large wavelength, and in general there may be non-adiabatic modes as 
well. These conclusions apply in all eras and whatever the constituents of 
the universe, under only a mild technical assumption about the wavelength 
dependence of the field equations for large wave length. In the absence of 
anisotropic inertia, the perturbations in the adiabatic modes are given for 
large wavelength by universal formulas in terms of the Robertson-Walker 
scale factor. We discuss an apparent discrepancy between these results and 
what appears to be a conservation law in all modes found for large wavelength 
in synchronous gauge: it turns out that, although equivalent, synchronous 
and Newtonian gauges suggest inequivalent assumptions about the behavior 
of the perturbations for large wavelength. 
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I. INTRODUCTION 



If observations are to be used to tell us something about inflation, then 

we need some way of connecting the properties of the cosmological fluctua- 
tions produced during inflation to the properties of fluctuations much closer 
to the present. Inconveniently, in inflationary cosmologies the era of inflation 
was followed by a period when the energy in scalar fields was converted into 
matter and radiation, and about this process we know essentially nothing. 
Subsequently there may have been other periods about which we are equally 
ignorant, such as the often-hypothesized era with temperatures between 10^^ 
GeV and 10^^ GeV, when supersymmetry may have become broken by un- 
known strong forces. These mysterious eras occurred when fiuctuations of 
cosmological interest were far outside the horizon, but this does not rule 
out some effect on the strength or even the wave-length dependence of these 
fluctuations.^ Therefore, in relating the cosmological fluctuations produced 
during inflation with those observed in the cosmic microwave background or 
in large-scale cosmic structures, it is essential to employ some sort of conser- 
vation law that is valid at large wavelengths whatever the details of cosmic 
evolution. 

In much work on fluctuations in cosmology, the conserved quantity is 
taken to be a quantity TZ related to the spatial curvature on co-moving spatial 
surfaces[l], given in Newtonian gauge by^ 

7e = + H5u . (1) 

The rate of change of TZ is given by a general formula[2] : 



7^ 




'H + 3HH\ ^ (H ^. 



(2) 



^By a mode being "beyond the horizon" we only mean that the physical wave number 
is much less than the expansion rate. This does not necessarily have anything to do with 
causality; indeed, the point of inflation is to make the true particle horizon radius much 
larger than the inverse expansion rate. 

^Here H = a/ais the expansion rate, with dots denoting ordinary time derivatives. In 
Newtonian gauge the perturbations to the gravitational field arc taken to be 5goo = —2$ 
and 5gij = —2a?"^5ij. Also 5p, dp, and 6u are the perturbations to the total energy density, 
pressure, and velocity potential in Newtonian gauge, while we use a bar to denote unper- 
turbed quantities like the unperturbed energy density p and pressure p. For simplicity we 
assume a vanishing unperturbed spatial curvature. 
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Here q is the co-moving wave number, 5a is a measure of the anisotropic 
stress,"^ and 

Thus TZ is conserved in the hmit of small wave number if and only if X = 
in this limit. 

The limit of small q is of some interest in itself, but its importance lies 
chiefly in the circumstance that those factors of q that arise from the fleld 
equations (rather than from the initial conditions) are always accompanied 
by factors of l/a{t), because it is only q/a{t) that is independent of the 
units chosen for the co-moving spatial coordinates x^. It is usually a good 
guess that terms in the perturbations proportional to such factors of q will be 
negligible if q/a{t)H{t) -C 1. Hence, although here we will study the behavior 
of the perturbations in the limit of small g, it is expected that this provides 
insight to the behavior of perturbations as a{t) increases. In particular, we 
expect that if X — > for g — > 0, and if the coefficient of q^ in the second 
term in Eq. (2) remains finite in this limit, then at any given epoch TZ will 
be small if q/aH is sufficiently small. In any mode for which TZ is non-zero 
in this hmit, the fractional rate of change of TZ will then also be small in this 
limit. 

Now, the quantity X vanishes in all modes and for all q when the per- 
turbed pressure p + 5p is a function only of the perturbed energy density 
p + Sp, as is the case in a universe dominated either by pure radiation or by 
pure cold matter, but not when both radiation and cold matter are impor- 
tant, and also not during inflation or in the curvaton model[4]. The quantity 
X does vanish for all modes in the limit 5 — > in the case of inflation with 
a single scalar fleld, but this is not true with several scalar flelds. Section 
II of this paper will show that in general, whatever the contents of the uni- 
verse, with only a mild technical assumption about the dependence on wave 
number of the fleld equations for cosmological perturbations in Newtonian 
gauge, these equations always have a physical solution for which X — > and 
TZ approaches a non-zero constant in the limit q ^ 0, though there may also 
be other modes for which TZ is not constant. In fact, there are always at least 

^The quantity Sa is defined by writing for scalar perturbations as gijp + didjSa. 
In this formalism, p includes effects of bulk viscosity, while iq6u is the velocity of energy 
transport, and so includes effects of heat conduction. For this formalism, see ref. [3]. 
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two solutions with X = and 71 constant in the hmit q ^ 0: one solution 
for which 71 0, and another with 71 = 0. These solutions will be illustrated 
in Section III for the case of inflation with any number of interacting scalar 
fields. The existence of such solutions is well known in special cases, but I 
do not know of a previous general proof of their existence. 

Solutions of this sort arc usually called adiabatic, even in contexts where 
thermodynamics has no relevance.'^ As we will also see in Section II, in 
theories in which the energy- momentum tensor is the sum of a number of 
tensors Tj^ for a set of fluids labeled /, we have the stronger result that for 
g = in the adiabatic modes, the perturbations in each fluid satisfy 

5pf Spf 5p 5p 

Pf Pf P P 

In the special case where the unperturbed energy-momentum tensors are 
separately conserved, we also have Pf — —3H{pf+pf), in which case Eq. (4) 
implies that the ratios 5pf/{pf +pf) are equal: 

5pf 6p 

Pf+Pf p + p' 

which is often taken to be what is meant by an adiabatic perturbation. 

Things appear very different in synchronous gauge. As shown in Section 
IV, when we take the limit of vanishing wave number in the field equations of 
synchronous gauge, we find that these equations respect the conservation of a 
quantity A in all modes whatever the contents of the universe, provided only 
that none of the perturbations in synchronous gauge blow up in the limit q — >■ 
0. At first sight, this presents an apparent paradox. All gauges are equivalent, 
so how can there be a quantity that is conserved for zero wave number in 
all modes during all eras in synchronous gauge and no equally universal 
conservation law for zero wave number in Newtonian gauge? We will find the 
answer to be that when we speak of the limit g = we mean different things 
in different gauges. Though mathematically equivalent, synchronous gauge 

^Sometimes the non-adiabatic solutions of the field equations are called isocurvature 
perturbations. This is a misnomer, because even for q = it is only possible for TZ to 
have the constant value zero if X = 0. Nevertheless, as wc will see in Section IV, in 
synchronous gauge there is a sense in which the solutions that do not correspond to the 
adiabatic solutions of Newtonian gauge can indeed be regarded as isocurvature modes. 
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(4) 
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and Newtonian gauge suggest different hypotheses about how perturbations 
behave in this hmit, leading to different conditions for the validity of the 
conservation law ior q — 0. 

In some work on cosmological fluctuations, a quantity related to the 
spatial curvature on spacelike surfaces of constant energy density, is used in 
place of TZ. It is defined in Newtonian gauge by [5] 

C = - H5p/p = + 5p/3{p + p) . (6) 

By taking a suitable linear combination of the time-time and time- space 
components of the Einstein field equation and the part of the space- space 
component proportional to 5ij, one can derive a general constraint [6], 

aHp - 3Ha' {p + p)5u+ v]/ ^ , (7) 

so that 

n-c, = { * (8) 

and therefore in all modes C, ^TZin the limit g ^ 0. 

II. ADIABATIC MODES IN NEWTONIAN GAUGE 

We consider a general cosmological model, based on the Einstein field 
equations, supplemented with whatever other equations are needed to give 

the energy density, pressure, and velocity potential perturbations in terms of 
independent dynamical variables, and the field equations satisfied by those 
variables. We will demonstrate the general existence of a pair of adiabatic 
solutions of these field equations in Newtonian gauge: one with TZ ^ and 
constant in the limit g — > 0, and the other with 7?. = in this limit. In order 
to draw these conclusions without specifying the formulas for dp, Sp, and Su, 
we use a trick, based on the fact that although there is no remaining gauge 
symmetry in Newtonian gauge for q ^ 0, for q = there is a remnant gauge 
symmetry of the field equations in Newtonian gauge, which makes it easy to 
find exact general solutions of the field equations for g = 0. Not all of these 
solutions are physical. To be physical, such a solution must be the limit as 
g — of a solution of the field equations for g 7^ in at least a neighborhood 
of g = 0. For this to be the case, it is necessary for the g = solution 
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to satisfy certain conditions imposed by the Einstein field equations, which 
limit the physical solutions to a linear combination of just two independent 
adiabatic modes. We will then make a mild technical assumption about 
the dependence of the field equations on q, which guarantees that these two 
solutions are the limits as g — of physical solutions for g 7^ 0. 

Whatever the constituents of the universe, for q = the field equations 
in Newtonian gauge for scalar (i.e., compressional) modes will always be 
invariant under the gauge transformation induced by a redefinition of the 
time coordinate 

t^t + e{t), (9) 
and a re-scaling of the space coordinate 

x'^x'{l-X), (10) 

with e(t) an arbitrary infinitesimal function of time and A an arbitrary in- 
finitesimal constant. For scalar modes, Sgio is proportional to qi, and the 
part of 6gij not proportional to 6ij is proportional to qiqj, so both auto- 
matically vanish for g — > 0, and we therefore do not need to impose any 
conditions on e(t) and A to remain in Newtonian gauge for q — 0. Eqs. (9) 
and (10) provide the most general spacetime transformations of purely scalar 
perturbations that preserve the condition q = 0. 

Under this gauge transformation, the q = perturbations undergo the 
transformation 

^ + He- X, $^$-e, (11) 

Sp ^ 6p — pe , 6p 6p — pe . (12) 

Likewise, in the case of separate fiuids with energy density and pressure pf 
and p/ or separate scalar fields </?/, 

^Pf ^Pf - Pf^ , hi hi - Pi^ , (13) 

or 

5(pf — > 5(fif — (ffC . (14) 

It follows that there is always a solution of the Newtonian gauge field equa- 
tions for g = 0, in which: 

^ = He-X, $ = -e , (15) 
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5p = —pe , Sp = —pe , 

and for several fluids 



(16) 



or several scalar fields 

5(pf = -(^_^e , (18) 

where e{t) is an arbitrary function of time and A is an arbitrary constant. (It 
is not necessary for this that the energy-momentum tensors of the individual 
fluids or scalar flelds be separately conserved; all we need is that they are 
tensors.) 

Of course, for general e{t) and A this is just a gauge mode. For it to 
have any physical signiflcance, it must satisfy certain conditions that allow 
it to be extended to the case of non-zero wave number. In particular, the 
part of the space-space component of the Einstein field equations that is not 
proportional to Sij takes the form (with Sa the anisotropic stress): 

qiqj {^-■qf) = -SnGqiqjSa , (19) 

so this equation disappears for q = 0. In order for the solution (15)-(18) of 
the field equations to be extended to q ^ 0, we must have 

$ = ^ - SnGSa , (20) 

and therefore e{t) and A must satisfy the condition 

e + He = X- 87rG6a , (21) 

which can always be satisfied by a suitable choice of e{t). 

There is another equation that disappears for q = 0: for scalar modes the 
space-time component of the Einstein field equations reads 

?i ( - 2^ - 2H^) = SttG (p + p) qM = -2Hqi5u (22) 

where 5u is the velocity potential in Newtonian gauge, which does not appear 
in the equations for g = 0. Hence in order for the solution we have found to 
be extended to g 0, this solution must also have a velocity potential given 

by 

HSu = ^ + , (23) 
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or, using the result (15), 

6u = e . (24) 

This agrees with what would be found from the gauge transformation induced 
by the coordinate transformation (9). Likewise, with several fluids, 

Suf = e (25) 

Prom Eqs. (16) and (24) it follows that Sp — 3H {p + p) 5u, so the constraint 
equation (7) is automatically satisfied for 5 = by this solution. By inserting 
Eqs. (15), (16), and (24) in equations (1) and (6), we now find that for g = 0, 
this solution has 

7^ = C = A , (26) 

so TZ and C, are indeed constant and equal for this solution in the limit g ^ 0, 
as was to be shown. They are also non-zero, as long as we take A 7^ 0. 

Now we have to ask what additional conditions are needed to ensure that 
this solution is the limit as g — > of a solution with g 7^ 0. In general, any 
closed set of linear homogeneous ordinary differential equations for a finite 
set of dependent variables can be put in the first-order form 

yn{t)+Y.Cnm{t)ym{t)^Q. (27) 
m 

(If some of the original set of field equations involve derivatives of higher 
than first order, we can still write the equations in the form (27) by including 
some derivatives of the field variables among the Umif)-) This has the general 
solution ^ ^ ^ ^ 

Vraito) (28) 



l/„(t)=^[T|exp(-£c(t')c?t')} 



nm 



with to arbitrary, and with T denoting a time-ordered product defined by a 
power series expansion of the exponential, which for finite matrices is always 
convergent. The initial conditions may be subject to constraints like Eq. (7), 
which can be written 

ECn?/n(io)=0. (29) 

n 

(Eq. (7) is such a constraint, because the equations of energy and momentum 
conservation and the gravitational field equation \I' + if$ = H5u imply that 
the left-hand side of Eq. (7) is time-independent.) 
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Our "mild technical" assumption is that, as long as the Einstein equa- 
tions (19) and (22) are written instead in the stronger form (20) and (23), 
the matrix elements Cnmit) and the constraint coefficients c„ are continuous 
functions of q in at least a neighborhood of g = 0. In this case the yn{io) 
that satisfy Eq. (29) and the matrix in Eq. (28) will also be continuous in 
g, so that any solution of Eq. (27) for g = can be extended to a solution 
for g 7^ in a neighborhood of g = by using Eq. (28) and (29) with the 
values of C{t) and c for g 7^ 0. The next section shows the validity of this 
assumption in one illustrative example: inflation with several scalar fields 
and an arbitrary potential. Because g generally enters the field equations 
and constraint equations in such a simple manner, we expect that this as- 
sumption will always be satisfied, and in any case it is easy to check for any 
specific model. For g 7^ 0, there is no remaining gauge freedom in Newtonian 
gauge (as there is for g = 0), so the adiabatic solution found in this way will 
be the limit as g — of a physical solution, not a mere gauge mode. 

We expect the anisotropic stress coefficient da for a wide class of theories 
to be some hnear combination of 5m, 5p, and so that for the solution 
(15)-(17), (24)-(25), 5a may be written os 5a — eE, with E(t) depending 
only on unperturbed quantities. (For instance, a non-zero shear viscosity 77 
in an imperfect fluid gives 5a = ~2rj 5u, so here E = —277.) In all such cases, 
Eq. (21) has the general solution 



and the lower limit on the integral in Eq. (30) is arbitrary. 

There is also a second mode, corresponding to the possibility of shifting 
the lower limit of the integral in Eq. (30), for which e(t) goes as e{t) oc l/a{t). 
Since shifting the lower bound on the integral in Eq. (30) has no effect on 
the value (26) of TZ and C, this solution has TZ — ( — 0. 

In the special case of vanishing anisotropic stress we have 5a = 0, so here 
$ = \I', and a{t) is just the Robertson- Walker scale factor a{t). The general 
solution of Eq. (21) is then 




(30) 



where 




(31) 




(32) 
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with an arbitrary lower limit. This eventually increases in absolute value as t 
for Robertson- Walker scale factors that grow as any power of t, while in the 
other mode e{t) oc l/a{t) decreases with time. Inserting the result (32) in 
Eqs. (15)-(18) gives explicit results for the perturbations in the gravitational 
field and various pressures and energy densities as functions of time. 

The results presented in this section can be interpreted in terms of what 
Liddle and Lyth in ref. [1] call a "separate universe" picture, which in one 
form or another has been used since the beginning of inflationary theory to 
deal with cosmological fluctuations in the case of a single scalar held. For 
instance, Bardeen, Steinhardt, and Turner in ref. [5] gave what they called a 
'heuristic argument' that in this case any portion of the universe that is larger 
than the horizon 1/H but smaller than the physical perturbation wavelength 
a/q would have to look like a separate unperturbed universe, with cp + 5(p 
following the unique evolutionary path of the scalar fleld, and with all of 
these separate universes therefore the same except for a variation in the time 
at which the scalar fleld satisfles some spcciflc initial condition a few Hubble 
times after horizon exit. As pointed out by Bardeen et al, it follows then 
that 5p/p — 5p/p, and hence X — Q, for q/a <^ H. 

There is a potential problem with this sort of argument, that there are 
two flelds involved, the inflaton and the gravitational fleld, so that different 
separate universes might have different ratios of these flelds. The argument 
of Bardeen et al. was formulated in a gauge in which it is unnecessary to 
consider fluctuations in the gravitational fleld, but it applies also to Newto- 
nian gauge, because in this gauge the constraint (7) allows the gravitational 
potential to be expressed in terms of fluctuations in the scalar fleld. But as 
we have seen in this section, in Newtonian gauge it is necessary not only to 
allow shifts in the time at which the scalar fleld reaches some given value 
after horizon crossing, but also to re-scale the co-moving coordinates used 
in each separate universe. In synchronous gauge there is no constraint like 
Eq. (7) that allows us to express the gravitational fleld in terms of the scalar 
fleld, and so, as we will see in Section IV, the solutions even for inflation with 
a single scalar fleld do not satisfy X = in the limit g — > 0. 

There is another potential problem, that the equation of motion of the 
scalar fleld is a seconc?-order differential equation, so that there are two in- 
dependent solutions whose relative coefficients may vary from one separate 
universe to another. Bardeen et al. and other authors avoid this problem by 
assuming that the scalar field experiences a period of "slow roll" infiation. 
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in which the differential equation satisfied by the scalar field is of first order, 
to a good approximation. We have not had to make this assumption, for a 
reason already pointed out by Guth and Pi [7]: the Wronskian of these two 
solutions decays rapidly after horizon crossing, so that it is as if there were 
only one independent solution. (Guth and Pi considered the case of H con- 
stant, but even with a time-dependent H the Wronskian still decays, though 
not precisely exponentially.) 

In any case, it has always been clear that such "separate universe" argu- 
ments do not rule out non-adiabatic solutions in the case of several scalar 
fields, where ratios of the scalar fields may vary from one "separate universe" 
to another. The results of this section may be interpreted as the statement 
that in this and all other cases it is always possible to find an adiabatic solu- 
tion of the field equations in Newtonian gauge in which the separate universes 
appear the same, except for a shift in the time coordinate and a re-scaling of 
the co-moving space coordinates. 

III. AN EXAMPLE: MULTIFIELD INFLATION 

For illustration, and to confirm the reasoning of the theorem of the pre- 
vious section in a case where X does not vanish for all modes, let us consider 
the case of inflation with an arbitrary number of scalar fields <^/, and with a 
general potential V that may include interactions among the various scalars. 
The energy-momentum tensor of the scalar fields has the perfect-fluid form, 
so here (7 = 0, and $ = ^f. The field equations in Newtonian gauge are 

if ^H-^ ^A-KGY,ip^5^f , (33) 
/ 

% + iHS^, + E g| + = -2*^ + . (34) 

and the constraint (7) is here 

(^i7+^j^ = 47rG^(-^^50/ + ^/^/) . (35) 

We can write Eqs. (33) and (34) in the form (27) by taking the y„ to run 
over ^ and all 0/ and 0/, in which case the constraint (35) is of the form 
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(29). Here obviously Cnm{t) and c„ are continuous in g in a neighborhood 
of g = 0; in fact, they are just hnear functions of q^. Hence any solution of 
Eqs. (33)-(35) that we find ior q — can be extended to a solution for q 0. 

Let us try for a solution for g = in which all of the individual velocity 
potentials —S(ff/(ff are equal, so that 

5(pf — —(ffSu , (36) 

with the common value satisfying 

5u^-^ . (37) 

Using the time-derivative of the unperturbed scalar field equation 

<^/ + 3/^<^/ + ^ = 0, (38) 

we can put Eq. (34) for g = in the form 

H5u + H5u + 5u^0. (39) 

Also, the gravitational field equation (33) now reads \E' + H'^ = H6u, which 
Eq. (39) guarantees is automatically satisfied by the ^' given by Eq. (37). 
The general solution is 

Su^-fadt, ^^HSu-X, (40) 
a J 

(with A an arbitrary constant), just as we found above in Eqs. (15), (24), 
and (32). The perturbations to the energy density and pressure of the /th 
field here are 

5pf = — + (p5Lp + -^^Sifif = — + Siijifij: — PfSu = —pfdu , (41) 
and 

5pf — —^ifi + (fSip — -^^5(pf — —y^ + Siijipf — pfdu — —pfSu , (42) 

so this mode is adiabatic, in the sense that X — > for q ^ 0. Inserting 
Eq. (40) in Eq. (1) gives again 71 = X. 
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Once again, because of the freedom to shift the lower hmit of the integral 
in Eq. (40), there are two adiabatic modes here, the second with 5u oc 1/a 
and 7^ = 0. 

For a single scalar field, Eqs. (33) and (34) are a third-order set of dif- 
ferential equations, and therefore have a third independent solution. The 
third solution can also be found explicitly, and turns out to have IZ oc 1/ a^H 
for g = 0, so this solution is not adiabatic. However, this third solution is 
eliminated by the constraint Eq. (35), which as we have seen in the previ- 
ous section is automatically satisfied by any adiabatic solution, but is not 
satisfied by the non- adiabatic solution of Eqs. (33) and (34). For scalar 
fields Eqs. (33) and (34) have 2N + 1 independent solutions, of which two are 
adiabatic, and one is eliminated by Eq. (35), leaving 2N — 2 non-adiabatic 
solutions. 

IV. SYNCHRONOUS GAUGE 

We now turn to synchronous gauge. With zero unperturbed spatial cur- 
vature, the perturbed metric has components 

5'ii(x, t) = a'^{t)6ij + /iij(x, t) , 5(00 = -1 , Qio = , (43) 

with hij a small perturbation. We now assume for simplicity that the per- 
turbed energy-momentum tensor takes the perfect-fluid form 

T^iu ^ pgiMu + {p + p)ui^Ui, . (44) 

The unperturbed quantities p and p depend only on time, and the unper- 
turbed velocity four- vector has components — 1, u'' — 0. The normaliza- 
tion condition u^u^ = —1 then requires that the velocity perturbation 5?i[f-* 
is purely spatial. (A superscript (S) is used to denote perturbed quantities 
in synchronous gauge.) We consider only compressional modes, for which 
Su\^^ = d5u^^^ / dx^ . Then the relevant fleld equations for a Fourier compo- 
nent with wave number q are [8] 

^ (a^V) = -47rGa^ + 3(5p(^)) , (45) 

and 

+ ZH + = - (p + p) - a-\Hu^^^) . (46) 
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Here ip is a field employed in recent work using synchronous gauge [9] 

There is also an Euler equation that will be needed later in this section: 

d 



dt 



(48) 



From equations (45) and (46) together with the relation A7rG(^p + pj — —H 
it follows that 



where 

Here is the proof: Eq. (46) can be written 



(49) 
(50) 



H 



and it follows immediately from Eq. (45) that 



d{a^HiP) 
di 



Eliminating Hip from these two equations gives 

d{a^Hilj) 



dt 



= -4TTGa'H{Sp^''^+3Sp^'^)+4iTGa' [Sp^'^ + 3H(Sp^'^ + Sp^'^)]+q'HSu^'^ 



or in other words 



d_ 

di 



The quantity in square brackets on the left is not invariant under the gauge 
transformations that preserve the condition (43) for synchronous gauge, so 
instead we work with the related gauge-invariant quantity (50), for which 
Eq. (49) follows immediately. 
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As long as the velocity potential remains finite in the limit g — > 0, Eq. (49) 
yields a conservation law 

A = for g = 0. (51) 

This is true for all modes in all cases, including inflation with several scalar 
fields and for the transition from radiation to matter dominance. The con- 
servation of A in the limit q = can also be derived by simply perturbing 
a(t), pit), and the curvature constant K in the Friedmann equation, which 
gives 5K = -2A/3. 

By taking suitable linear combinations of solutions, it is always possible 
to arrange that for g = just one of a complete set has A ^ 0, while all 
the other solutions have A = 0. Examples are given in an appendix to 
this paper. Because of the connection of A with the spatial curvature, it 
is legitimate to call the solutions with A = isocurvature modes. When 
q is small but non-zero the isocurvature solutions usually have both A and 
A of order g^, so that Eq. (49) does not keep A for these solutions from 
undergoing large fractional changes. This does not vitiate the usefulness of 
the conservation law for initial conditions that give a physical perturbation 
in which all solutions make contributions with comparable coefficients. In 
this case, the contribution of the isocurvature modes to A may be rapidly 
varying, but at any given time they will be small as long as q is sufficiently 
small. The physical solution will have a rapid fractional variation in A only 
if the coefficient of the mode with A ^ for g = is suppressed, or if the 
coefficients of the isocurvature modes are enhanced. 

There is a simple relation between the quantity A introduced in this sec- 
tion and the more familiar quantity TZ discussed in Sections I-III. Given 
perturbations ^, Sp and Su in Newtonian gauge, we can find the pertur- 
bations 'tp, Sp^^\ and du^^^ in synchronous gauge from the transformation 
equations: 

iP = -3^ - 3^ (He) + (q/a) % , Sp^^^ ^Sp-ep Su^^^ ^Su + e, (52) 
where 

e = ^ . (53) 

(The possibility of shifting e by a constant term corresponds to the possibility 
of making gauge transformations that preserve the conditions for synchronous 
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gauge.) By applying these equations to the quantity (50), it is elementary 
to show that A is related to the quantity 71 defined in Eq. (1) by 



A = -q^TZ . (54) 

Thus for any finite q the fractional rate of change in TZ will be the same as the 
fractional rate of change in A. In some treatments of multi- field infiation[10] 
and in discussions of the curvaton model[4], it is simply assumed that the 
mode with 7^ 7^ and hence A ^ is somehow suppressed, which is enough 
to explain why these authors find a significant fractional change in TZ. But 
why more generally does the condition X — play an important role in 
establishing the conservation of 7?. for — > in Newtonian gauge, while 
there seems to be no similar condition needed for the conservation of A in 
synchronous gauge? 

As a first step toward resolving this apparent paradox, we note from 
Eq. (54) that that the limit as 5 —> of the perturbed quantities Sp^^^ and 
ip in synchronous gauge in the mode for which A ^ in this limit is not 
obtained by applying a gauge transformation to the perturbed quantities in 
the corresponding mode in Newtonian gauge for q = 0, since that would give 
A = for q = 0. We can go further, and show in general that for g = 0, 
the synchronous gauge solution corresponding to any adiabatic solution of 
the Newtonian gauge field equations (normalized to not diverge as g — > Oj 
has vanishing values not only for A, but also (up to a choice of a particular 
synchronous gauge) fort/j and the total density fluctuation 5p^^^ and velocity 
potential 5u^^\ 

The reasoning here is essentially the reverse of that used to prove the 
theorem of Section II. We will use the space-time component of the Einstein 
field equations in Newtonian gauge 

^ + H-^ = -47rG(p + p)Su = HSu . (55) 

We work in the limit 5 = 0, assuming that the sohition is normalized so that 
in Newtonian gauge all fluctuations remain flnitc in this limit. (As we shall 
see, this assumption is less innocent than it may seem.) Then for modes that 
for = are adiabatic in the sense that X = 0, Eqs. (1) and (2) give 

* = I (HSu) (56) 
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for q = 0. Combining this with Eq. (55) gives 



^ ^-Su. (57) 

Thus according to Eq. (53) we can adopt a particular synchronous gauge 
such that the transformation parameter e in Eq. (52) is 

e = -6u . (58) 

Using Eqs. (56) and (58) in Eq. (52) shows immediately that, for q — 0, 

ijj^O. (59) 

Furthermore, Eq. (55) together with the Newtonian gauge Euler equation 
supplies a general constraint equivalent to Eq. (7) for q — 0: 

-A7rG6p = 3HH6u (60) 

Eqs. (52) and (60) give the synchronous gauge density fluctuation 

-AnGSp^^^ = -A7rG[5p-ep] 

= 3HH5u + [-A7rG]5u [-3H{p + p)] 

= 0. (61) 

Finally, the velocity potential in this synchronous gauge is 

^Su + e^O. (62) 

Thus no synchronous gauge perturbation with non- vanishing values of ip 
or Sp^^^ or Su^^^ (apart from those that can be eliminated by a transformation 
to a different synchronous gauge), such as modes 1,2, and 3 of the radiation 
plus cold dark matter model of the appendix, can be the gauge transforma- 
tion of one of the q = adiabatic Newtonian gauge solutions. Rather, the 
synchronous gauge solutions ior q — with A a non-zero constant must be 
the gauge transformations of the terms of order in the adiabatic Newtonian 
gauge solution with TZ ^ 0, re-normalized by dividing by a factor g^.^ With 

^This is why it is possible for the quantity X not to vanish in any mode for q = in 
synchronous gauge, as we find in the appendix in the case of inflation, while there are two 
modes in Newtonian gauge in which X — > for g — > 0, despite the fact that X is gauge 
invariant. It is not that X is different in the two gauges, but rather that the limit q 
means different things in synchronous and Newtonian gauge. 
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this re-normalization of the synchronous gauge modes, as in the appendix, 
the conserved quantity A is not necessarily of order q"^, as would be expected 
from Eq. (54), but can have a finite limit for g — > 0, as we will find it does 
in the appendix. 

Now at last we come to the point. Working in Newtonian gauge, it is 
most natural to assume that, with an over-all normalization factor chosen 
so that TZ is finite and non-zero in the limit q = 0, all density fluctuations 
and velocity potentials as well as ^ are non-zero in this limit. Under this 
assumption, if the contribution of non-adiabatic modes is comparable to that 
of the adiabatic modes, TZ will undergo significant changes with time. Trans- 
forming this sort of solution to synchronous gauge, we have found above that 
the density fluctuations and the total velocity potential receive contributions 
of order (relative to the Newtonian gauge perturbations) from the adia- 
batic modes but of order unity from the non-adiabatic modes, so that A is of 
order g^, while A is also of order g^, and so A does suffer signiflcant changes 
with time. Or we can re-normalize the synchronous gauge fluctuations by 
an over- all factor of order 1/g^, in which case A and the density fluctua- 
tions and velocity potentials receive contributions of order unity for g = 
from the adiabatic modes, as in the appendix, while the contribution of the 
non-adiabatic modes to the total velocity potential if present is enhanced by 
a peculiar looking factor of 1/g^, giving both A and A non-zero limits for 
g ^ 0. 

On the other hand, working in synchronous gauge, it is most natural 
to assume that, with an over-all normalization factor chosen so that A is 
flnite and non-zero in the limit g = 0, all density fluctuations and velocity 
potentials as well as ip are flnite in this limit. Under this assumption, it makes 
no difference whether the contribution of non-adiabatic modes is comparable 
to that of the adiabatic modes; even if it is, A will undergo no significant 
changes with time. Transforming this sort of solution to Newtonian gauge, 
one finds that the density fiuctuations and the total velocity potential receive 
contributions of order 1 /g^ from the adiabatic modes and of order unity from 
the non-adiabatic modes, so TZ is of order 1/g^ while its rate of change is only 
of order unity. Or we can re- normalize the Newtonian gauge fluctuations by 
an over- all factor of order g^, in which case TZ and the density fluctuations and 
velocity potentials receive contributions of order unity for g = 0, while the 
contribution of any non-adiabatic modes is suppressed by a peculiar looking 
factor of g^, giving TZ a zero rate of change for g = 0. 
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So which is right? The issue is not the over-all normalization of the total 
perturbations, but the relative magnitude of its adiabatic and non-adiabatic 
terms in the limit g — > 0. There is nothing about either gauge that makes it 
a more reliable guide to our intuition about this than the other. 

It is generally expected that inflation with several scalar fields the gen- 
eral solution docs not have TZ approaching a constant for increasing a{t), in 
agreement with what would be expected from the behavior for g — > sug- 
gested by Newtonian gauge but not synchronous gauge. But there are cases 
of multi-field inflation in which A and hence 71 do approach constants as a{t) 
increases, as would be expected from the behavior for g suggested by 
synchronous gauge but not Newtonian gauge. One case is a potential given 
by a sum of exponentials[ll]. 

V ^J29n exp(- A„(/7„) (63) 

n 

Another is a potential of the form 

(64) 

with F an arbitrary function. It would be interesting to characterize the 
general class of potentials for multi-field inflation for which A and TZ approach 
constants as a{t) increases. 

APPENDIX: LONG- WAVELENGTH SOLUTIONS IN 
SYNCHRONOUS GAUGE 

In this appendix we will study several examples of calculations for zero 
wave number in synchronous gauge, to exhibit both solutions with A ^ 
and those with A = 0. All quantities here will be in synchronous gauge, so 
we will drop the label (S). 

As a first example, consider inflation with just a single real scalar field 
(p — (f{t) + Sip{x, t), and potential V{(p). As is well known, the unperturbed 
pressure and energy density are 

p^^^' + V{cp), p=l^^-V(^), (65) 
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from which we find the equation of motion of the unperturbed scalar field 



(p + 3//0 + V'{^) = . (66) 

The perturbations to the energy density and pressure are 

5p — (p5(p + V'{0)5ip , 5p — (p5(p — V'{(f)5(f . (67) 

Also, the perturbed velocity potential is 

5u — —5(p/(p . (68) 

The field equations (45) and (46) for the Fourier component of the pertur- 
bations with wave number q here take the form 

A (a^^^ ^ -AnGa' {^0 Sep - 2V'{<f) Sep) , (69) 

(JjV 

S(fi + 3H6(p + V"{(fi) 6(p + a-\^S(p = V , (70) 



where 



a ^^ 



Vm . (71) 



The Euler equation (48) gives no new information here. 

There is a gauge mode, with ip = r^p and ip = t{3H — q'^/a^), where r 
is an arbitrary time-indcpcndcnt function of q. Knowing this solution allows 
us to reduce the degree of equations (67) and (68) from three to two, in 
agreement with the number of physical solutions found in Newtonian gauge 
in Section III. We introduce time-dependent functions / and g by writing 

dip = f0, ^={f + g) {SH - q'/a') . (72) 

Equations (67) and (68) then become a second-order set of equations for the 
gauge-invariant quantities / and g: 

f + 3Hf + {H/H)f = -{3H- qya')g , (73) 
{SH - q^/a^)g + (q>HH + 3H)g=[H + q^/a^)f , (74) 
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in which the gauge mode appears in the possibihty of adding a constant to /. 
These equations can be solved exactly for q = and H{t) arbitrary. There 
are two physical solutions: 



Mode 1: 



^(1) U) = —. /* a(t')dt' , (t) = -. 

Mode 2: 



Hit) , , „ , 
1 ^ / a{t')dt 



a{t) Jo 



t 

'\ jj-i 



(75) 



^(2) (^) = L_ , /(2) ^t) = . (76) 

(The lower limit on the integral over t' is arbitrary; changing it just amounts 
to adding some of mode 2 to mode 1.) 

Equation (50) gives the values of ^4 for g = in these two modes as the 
constants 

= 1 , ^2 = , (77) 

even though neither of these solutions satisfies the adiabatic condition X = 0. 
A general mixture of modes with coefficients Ci and C2 will have A — ci for 
q — 0, provided C2 does not blow up in this hmit. With this proviso, the 
conservation of A allows the value of Ci that is calculated for a given inflaton 
potential to be used to find the strength of the non-isocurvature mode at later 
times. However, if the value of C2 for the physical solution found after horizon 
crossing went as c/q'^ for g — > 0, while Ci remained finite, then Eqs. (49), (66), 
(70), and (74) would give 

A ^ C2i2 ^ ^^I^ for q^O, 
a\t)H{t) ^ 

and there would be no useful conservation law even for g = 0. As discussed 
in Section IV, this is just what we would expect if we assumed that, with 
an over-all normalization factor chosen so that TZ is finite and non-zero in 
the limit g = 0, all density fluctuations and velocity potentials in Newtonian 
gauge as well as ^ are finite and non-zero in this limit. 

For another example, we consider a later epoch, when the dominant con- 
stituents of the universe were radiation and cold dark matter. (For simplicity. 
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we are neglecting the baryon density compared with the density of cold dark 
matter, but supposing that there are still enough baryons to keep the radia- 
tion in thermal equilibrium, and we are ignoring the effects of free-streaming 
neutrinos.) We adopt a particular synchronous gauge in which the cold dark 
matter is at rest. The field equations then are Eqs. (46) and (48) for the ra- 
diation energy density perturbation 6pf> and velocity potential 5uji\ Eq. (46) 
for the cold dark matter density p^] and Eq. (45) for ■i/', with the total energy 
density and pressure appearing on the right hand side: 



5 PR + mSpR = -(4/3)pr(V^ - a-^q^6uB) , 



(78) 



4^ 
dt 



a^PRSuR = -a^Sp 



R 



5pD + SHSpD = -poi^ , 



d_ 
dt 



{a' 



-iTTGa^ {25pR + 6pd) 



(79) 
(80) 
(81) 



The unperturbed radiation and dark matter densities go as and a~^, 
respectively. It is convenient here to normalize a so that a = 1 when pR — po, 
so that 

PR = PEQa~'^ , Pd = PEQa~^ , (82) 

where peq is constant. 

Equations (78)-(81) are a fourth-order system of differential equations, 
so there are four modes, all of which are physical because the gauge has been 
fixed by choosing Sud — 0. For q — 0, they take the form 

Mode 1: 



Su 



(1) 

R 




16 + 8a - 2a^ + a^J , 



f 



4pEQ V S^ttGpeq J vT+a 
3 

16 + 8a - 2a^ + 
VTTa ' 




32 + 8a - a-" 



(83) 
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Mode 2: 



Mode 3: 



Mode 4: 



1 

a 



7 



(2) 



4p£;Q V 8^GPEQ J vTTa 
3 



SnGpEQ 



4 + 3a 



.(3) 



5u^^ oc a . 



6u 



(4) 
R 

(4) 



a 



+ 4a - - 8\/TT 



7 



da 



^Pd 



V ^T^GpEQ J a/1 + a 

3 



1^8 + 4a - 8Vl + aj , 

|^(4 + a)A/rTa-4-3a^ 



(84) 



(85) 



(86) 



The lower bound on the integrals in the formulas for 6uji in modes 1,2, and 4 
are arbitrary; changing this lower limit in any of these integrals just amounts 
to adding some of mode 3 to that mode. 

Note that modes 1,2, and (trivially) 3 are adiabatic, in the sense that 

SPD Spn ^g^^ 



PD 



PR + PR ' 



(and so X = 0) while mode 4 is not adiabatic in this sense. The values of A 
for the four modes are 



Ao = A:> = Aa = 0. 



(88) 
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Thus modes 2 and 3 are both adiabatic and isocurvature. An arbitrary 
mixture of modes will have A constant unless the coefficients of modes 2, 3, 
or 4 blow up as l/gr^ hmit ? — > 0, which will be the case if the fluctuations in 
the non-adiabatic modes in Newtonian gauge have non-zero limits for g — > 0. 
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